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Abstract. The restricted version of the Hilbert 16th problem for qua- 
dratic vector fields requires an upper estimate of the number of limit 
cycles through a vector parameter that characterizes the vector fields 
considered and the limit cycles to be counted. In this paper we give an 
upper estimate of the number of limit cycles of quadratic vector fields 
"(T-distant from centers and K-distant from singular quadratic vector 
fields" provided that the limit cycles are "5-distant from singular points 
and infinity". 

1. Introduction and statement of the main result 

Hilbert 16th problem asks (see [6]): what may be said about the number 
and location of limit cycles of a polynomial vector field of degree n in the real 
plane ? The main contributions in this direction were the works of Ecalle [5] 
and Ilyashenko [7] who proved that any polynomial vector field has finitely 
many limit cycles, and also the work of Llibre and Rodriguez [13] who showed 
that any finite location of limit cycles is realized by a polynomial vector field 
of a convenient degree. But the complete answer to Hilbert 16th problem 
is unknown even for n = 2. Even the existence of an uniform upper bound 
of the number of limit cycles for quadratic vector fields (polynomial vector 
fields of degree 2) is not yet proved. Limit cycles of a quadratic vector field 
may surround only one singular point, and it is of type focus (for more 
details see [H]). Moreover, quadratic vector fields have at most two foci 
(see again |14|). Limit cycles surrounding the same singular point form a 
nest. Recently Zhang Pingguang [T51 [TBI E] proved that only one nest of 
quadratic vector field may have more than one limit cycle. 

The restricted version of the Hilbert 16th problem for quadratic vector 
fields allows us to introduce a vector parameter that characterizes the vector 
field and the limit cycles to be counted. The upper bound for the number 
of limit cycles is expressed through this parameter. 
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In this paper we give an upper estimate of the number of limit cycles of 
quadratic vector fields "a-distant from centers and K-distant from singular 
quadratic vector fields" provided that the cycles are "<5 -distant from singular 
points and infinity". The precise sense of assumptions in quotation marks 
is explained below. The upper estimate mentioned above depends on a, k 
and 5. 

1.1. Normalized quadratic vector fields. We consider quadratic vector 
fields with a focus point and estimate the number of limit cycles that 
surround this point. The system has the form 

(1) z = [iz + Az 2 + Bzz + Cz 2 , 

where [i,A,B,C are complex coefficients. Rescaling: z i— > cz and t c't, 
c£C,c'eR brings it to 

c 2 

z = c'Uz + Acz 2 + Bczz + C—z 2 ). 

c 

Hence, after an appropriate normalization, we can take in (TjQ): fJ, = X% + 
i, max(|^4|, |C|) < 1, \B\ < 2. Moreover, the normalized tuple (A,B,C) has 
the form: either A = 1 and \B\ < 2, |C| < 1, or B = 2 and \A\ < 1, \C\ < 1, 
or C = 1 and \A\ < 1,\B\ < 2. The reason for distinguishing B will be seen 
later. 

To summarize, the normalized quadratic vector field has the form: 

(2) z = fj,z + z 2 + Bzz + Cz 2 , \B\ < 2, |C| < 1, 
or 

(3) z = fiz + Az 2 + 2zz + Cz 2 , \A\ < 1, |C| < 1, 
or 

(4) z = fiz + Az 2 + Bzz + z 2 , \A\ < 1, \B\ < 2, 
with jj, = \\ + i. 

Moreover, we consider that Ai > 0. If not, we reverse the time and make 
a symmetry z ^ z. 

The tuple of parameters (Ai, A, B, C) with {A, B, C) normalized as above 
is denoted by A and the corresponding vector field (and equation) is denoted 
by v\. The space of all these A's is denoted by A. It is homeomorphic to 
the glued union of three copies of M + x D 2 x D 2 , where M + = [0,oo) and 
B 2 = {z G C : \z\ < 1}; the gluing maps identify the boundaries of the cells 
R + x I 2 x D 2 ; we will not need these maps below. 

1.2. Center conditions. Center conditions for quadratic vector fields are 
well known; see the works of Dulac [4], Kapteyn |11} I12j. Bautin [1]. In the 
complex form they were found by Zoladek [18], see the next theorem. We 
will use the latter form of the center conditions. 
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Theorem 1. A quadratic vector field (JXJ) has a center at zero if and only if 
the following holds: 

= Ai = 0, 
= Im(AB) = 0, 

= lm[(2A + B)(A- 2B)BC)_= 0, 
= Im[(2A + B){\B\ 2 - \C\ 2 )B 2 C] = 0. 

Definition 2. A normalized quadratic vector field is called a-distant from 
centers provided that 

4 

(6) £bi(A)|>a. 

The set of all A G A /or which v\ is a-distant from centers is denoted by 
A(o-). 

1.3. (5 tame limit cycles. Now for any 5 > we define 5-tame limit cycles 
of a normalized quadratic vector field. Note that the normalization of a 
quadratic vector field provides a scale on the phase plane. Thus the following 
definition makes sense. For any 5 E (0, 1) and any A 6 A denote by B(X,5) 
the disc \z\ < <5 -1 in C minus all the open ^-neighborhoods of the singular 
points of v\, both real and complex, except for the point 0. 

Definition 3. A limit cycle of a normalized vector field is called <5-tame 
provided that it belongs to B(X,S). 

1.4. Singular quadratic vector fields. A quadratic vector field with a 
focus at the zero and a line of singular points not passing through zero is 
called singular. Such a quadratic vector field can be written as 

i = fizl(z), 

where l(z) is a real polynomial of degree 1 of the form l(z) = az + az + 1. 
After normalization, this equation has the form 

(7) i = [iz + z 2 + ^zz := vJz), 

where = X± + i. The s of v s is for a singular quadratic vector field. Any 
normalized quadratic vector field close to a singular one has the form 

(8) v = v s + u, u = bzz + cz 2 ; 

here v s is singular, b and c are small. In this expression for v, its coefficient 
B may be greater in modulus than 1 but smaller than 2 because b is small. 
Still the normal form of v is ( }2"]) . That is why B is distinguished in the 
definition of the normal form. To find decomposition (|8|) for a vector field v 
in the form ([2]), take v s as in ([7]) and u as in ([8]) with coefficients: 

b = B-^, c = C. 



(5) 
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Let || • H2 denote the Li norm on a circle. Then 

||r- 2 u||| = |6| 2 + |c| 2 . 

Definition 4. A quadratic vector field is K-distant from the set of singular 
quadratic vector fields if \ \r~ 2 u\\2 > k in (|8|). 

1.5. Main result. 

Theorem 5 (Main Theorem). For any {5, a, k} C (0,0.1), the number of 
5-tame limit cycles of a normalized quadratic vector field which is a -distant 
from centers and k -distant from singular quadratic vector fields is no greater 
than 

H(2,5,a,n) = |loga|exp(exp(10 25 ^ 31 ^ 2 )). 

This estimate is irrealistic but this is the only known estimate of this 
kind. 

This paper is the first in a series of papers aimed to estimate the number 
of 5-tame limit cycles of quadratic vector fields. In a subsequent paper we 
prove that for k sufficiently small: n < kq(5,o-), the vector field ([8]) has only 
one 5-tame limit cycle. A similar result, without a quantitative estimate on 
the value of Ko(5,a), is obtained in the preprint [3]. 

1.6. Growth— and— Zeros Theorem. Limit cycles correspond to the fixed 
points of the Poincare map . For normalized quadratic vector field v\ con- 
sider the Poincare map P\ of a segment of a positive semiaxes M + with the 
left endpoint into M + ; the right endpoint will be specified later. 

The number of the fixed points of this map will be estimated with the 
use of the theorem named in the title of the subsection; for its proof see [7J, 

ma. 



Theorem 6. Let U C C be a connected and simply connected domain and 
K C U be a path connected compact set. Let D be the internal diameter of 
K, and 

gap (K,U) :=p(K,dU)>e, 
where p(K,dU) = min^Ja — b\. Let f : 0U — > C be a holomorphic func- 
tion. Then 

(9) #{> £ K I f(z) = 0} < B KiU (f) exp {^-\ , 

where 

(!0) B K,u{f) = log rjr • 

max K I /| 

As usual 0U denotes the closure of U. The expression Bx,u{f) is called 
the Bernstein index of f for K, U. The exponential in (J9j) is called the 
geometric factor. We will often write: 

M = max l/l, m = max 1/1. 
U K 
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This theorem will be applied to bound the number of zeros of the displace- 
ment function 

fx = P\~ id 

of the Poincare map P\ of v x \ these zeros correspond to limit cycles of v x . 
There are the following steps in the application of this theorem: 

- choice of K and finding the lower estimate for m = max \f\\\ 

- choice of U and finding the upper estimate for M = max \f\\. 

2. The lower estimate of the maximum of the displacement 

2.1. Normalized quadratic vector fields in polar coordinates. To 

write equation ([I]) in polar coordinates (r, 9) note that 

n ogz y = t +i e = ^. 
r z 

Hence, 

r = rRe ^ = r(Ai + rf x (9)), 

t 11 ) ■ viz) 

9 = lm ^ = l + rg x (9), 

where f x and g x are trigonometric polynomials of degree 3: 
h(6) = Re h x (9), g x {9) = Im h x (9), 

(12) h x (9) = Ae i9 + Be~ i9 + Ce~ 3i9 . 
For the normalized equations, |/t,\(0)| < 4. Hence, 

(13) \f X (9)\ <4, \g x (9)\ <4. 

2.2. Compactification. 

Lemma 7. If a system v x has at least one 5 -tame limit cycle, then \X\\ < 
4/5. 

Proof. Let Ai > f Recall that r < in B(A, <5). Then in B(X, 5), f > by 
(jlip and (|13p . Hence, the vector field has no limit cycles in B(\,5). □ 

2.3. Complex extension of the Poincare map near zero. We will com- 
plexify nonautonomous equation corresponding to the system (jlip making 
r complex and denoting it by w and keeping 9 real. We get: 

riA\ dw ^i+wf x {9) . . ! 

(14) — =w-— — :=F x {w,9), w£C, 9e§. 

d9 l + wg x (9) 

Recall that ||/a|| < 4 and \ \g x \\ < 4. When the norm is not specified, it is 
the C-norm of a function on the circle. 

For any value of Ai , we will find R and e in such a way that the orbit that 
starts in a cross-section D £ := {\w\ < e} x {0} keeps inside W := {\w\ < 
R} x 8 1 when 9 ranges over [0, 2ir]. We call this property (*) of (fT4l) . 
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Lemma 8. Equation (fTSjl satisfies property (*) for R = 0.01 and 

' 0.001 /or Ai G [0,0.1] 



(15) £ = 2e(A) 



Proof. The proof is based on the Gronwall inequality that measures the 
divergence of two solutions of a differential equation. To apply the classical 
Gronwall inequality to a differential equation with the complex phase space, 
we simply take the realification of this space. In case when one of the 
solutions is identically zero, the Gronwall inequality measures the norm of 
the other solution. For equation (|14p this inequality has the folowing form. 
Let 

L = max — — 

w dw 

and |w(0)| < e. Then the Gronwall inequality claims that 

(16) \w(0)\ < ee Le for € [0,2vr], 
provided that 

(17) se 27TL < R. 
To get an upper bound for L, note that 

flg v dF\ = Ai + 2wf x _ wjXi + ™/a) 

' dw 1 + wgx (l + wg x ) 2 9X ' 

Note that ||/ A || < 4, \\g x \ \ < 4. Hence, 

'0.2 for Ai < 0.1, 
^2Ai for Ai > 0.1. 

Now, (dU) yields Lemma EJ □ 



(19) L < 



2.4. Description of K x . Let V be the positive x semiaxis. Assume that 
system v x has no <5-tame limit cycles around the origin. Then Theorem Q] 
holds for this system. In what follows, we consider the opposite case. Let 
a(A) be the intersection point of the outmost tame limit cycle surrounding 
the origin with T. Let s x be the segment [0, a(A)], and e(A) be the same as 
in (El). 



Lemma 9 (First Main Lemma). For the set 

(20) K x = s x U D £[x) 
the following lower estimates hold: 

(21) m(A) := max \P\{w) -w\> 10' 2e a for Aj < 0.1 and 

(22) m(A) > 10^ 26/5 for Ai > 0.1. 

Note that these estimates do not depend on k. The lemma is proved in 
the next five subsections. 
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2.5. Proof of Lemma [9] for the case of a strong focus. In this sub- 
section when we say that the normalized quadratic vector field has a strong 
focus we mean that Ai > 0.1. 

To prove Lemma [9] in this case, we use the reversed Cauchy inequality for 
the first derivative: if / is holomorphic in a disc D e = {\w\ < e} x {0} and 
continuous on the boundary of this disc, then 

(23) max|/|>£|/'(0)|. 

For / = P\(w) — w, and in the case Ai > 0.1, we have: 

/'(0) = e 27rAl - 1 > 0.3e 27rAl 
By LemmaEl / is well defined in D e for e = 0.005e~ 4Al?r . Hence 
m > max l/l > 0.0015e _2Al7r > e~ 26/5 , 

D e _ 

where the last inequality follows from Ai < | and 6 < 0.1. This yields (|22p 
and proves Lemma [9] for Ai > 0.1. To prove this lemma for Ai < 0.1, we 
need first to study the case Ai = and then to perturb it. 



2.6. Seven jet of the Poincare map for linear part a center. The 

Poincare map for the point zero of the normalized quadratic vector field v\ 
may be decomposed in a Taylor series 

(24) F A H = ^a i (AV'. 

This series converges at least in a neighborhood of the form D° = {\w\ < ro} 
for a convenient ro > 0. Consider the case Ai = 0. For such A, the coefficients 
Oj(A) become functions only of (A,B,C) not necessarily normalized. 

Lemma 10. Let Ai = 0. Then for the decomposition (I24p . 
ai = 1, a 2 = 0, a 3 = a g 2 , a 4 = aig 2 , 

«5 = Po93 + 0192, a e = (3 2 93 + 0392, a 7 = 7054 + 7153 + 7252, 

where g 2 ,93, 54 are the polynomials from the center conditions (|5j) , ay , f3j , 7j 
are polynomials in the variables A,B,C, and Qo,/?o ; 7o are constant. More- 
over, on the set of A = (0, A, B,C) with the tuples A,B,C normalized we 
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have: 

M < 2, \g 3 \ < 30, |# 4 | < 36; 



|«o| = 27T, 

2-7T 2tt 

\Po\ = y > IAI ^ y (284 + 1087r) := Bl < 500 ' 

l7o| = N ^ ^( 5816 + 15367r ) := Cl < 500 ' 

vr (5019144 + 2565120vr + 345600tt 2 ) r , 

M < ^ ^ := C2 £ [4 • 10 4 , 10 5 ]. 

Expressions for the a 's, (3 's and 7 's are given in the appendix. 

Lemma flOl has been proved using the algebraic manipulator mathematica 
and the three normal forms for the quadratic vector fields. The algorithm 
is sketched in the appendix. 

2.7. Lower estimate: case of a linear part a center. Denote the nor- 
malized tuple A with Ai = by A' := (0,A,B,C). Recall that in ([20]) . 
e(A') = 0.0005. Let m(A) be the same as in ([2T]) . Recall that A(cr) appears 
in Definition [2j The next lemma is one of the main steps in the proof of 
Theorem [TJ 

Lemma 11. For the normalized A £ A(<r) with Ai = 0, we have: 

m(A) > 2 • 10 _23 cr := m . 
Proof. Let f\ = P\ — id. By Lemma CEl for Ai = 0, A = A', we have: 

/a(0) = f x (0) = /a(0) = 0. 
For vector fields cr-distant from centers, we will prove a lower estimate: 

K(A)| > mj(a), 

with nij explicitly written for at least one j G {3; 5; 7}. By Lemma [8] the 
function f\ is holomorphic in the disc \w\ < 0.001 := 2e(A'). Hence, there 
exists j E {3, 5, 7} such that 

m(A) > max |/a| > mj(a) ■ e(X') J . 

D e{X') 

The lower bounds for Oj are found in the following way. For a, (3 £ (0, 1) 
chosen later, the compact set Ao(cr) = A(cr) n {Ai = 0} is split into three 
parts S 2 , S3, S4 where 

S 2 = (\g 2 \ > aa), S 3 = (\g 2 \ + M > 0a) \ S 2 , S 4 = A (a) \ (S 2 U S 3 ). 

On Sj, |a2j-i(A)| is estimated from below. By Lemma [TOl on S 2 , \as\ > 
ctoaa. Let B\,C\, C 2 be the same as in Lemma [10J On S3 we have: 

05 = A)03 +/?102- 
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Hence, 



\a>5 |s 3 1 > Po(P - a)a - B ± aa = O (@ ~ a ^1 + jfjj °- 



If we choose a so small in comparison with j3 that 

BA < P 
do) ~ 2' 



(25) «(1 + ^<I 



then 

On £4 we have: 

As C2 > C%, we have: 



07 = 7o54 + 7i53 +7292- 



K k I > 7o(l - 0)<? ~ C 2 (3a = 7o I 1 - /3 I 1 + — ) ) a. 



If /3 is so small that 

(26) (,( 1 + £l)<l 



then 



Now, 



To 



70 7 " 2' 
Ms 4 |>— • 

\m |s 4 I > mm |a7|r > e(A J := m^a 

s 4 2 

\Tn I s 3 I > mm l a 5| r o — — o — ) := m 3 a 
S3 2 

I 711 - |s 2 I — mrn l a 3l r o — aoacre(A') 3 := m2<J. 

Due to Lemma OH inequalities ([25]), ([26]) hold for /? = 10" 5 ,a = 2 • 10" 8 . 
Again by Lemma [TO] Tcii > m 3 > m 4 > 2 • 10~ 23 . This proves Lemma [Til 

□ 

2.8. Proof of the First Main Lemma in case of the moderate fo- 
cus. Recall that mo is the lower estimate of maxo E \P — id\ mentioned in 
Lemma [11] Here we consider the case Ai G [mo, 0.1]. In this case, by 
Lemma [8] the displacement fx of the Poincare map is holomorphic in a disc 
| w | < £ = 0.0005. We have: 

1^(0)1 > e 2-no _ J > 27rmo . 

Hence, 

max |/a| > 0.003m,o. 

\w\<e 

This proves the First Main Lemma in the case considered. 
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2.9. Proof of the First Main Lemma in case of the slow focus. We 

consider here the last remaining case Ai G (0,mo], where mo is the same 
as in Lemma [TTj i.e. thq = 2 • l(T 23 o-. This case is treated small 
perturbation of the case Ai = 0. Consider two systems (fl4l) corresponding 
to Ai = and Ai € (0, mo] fixed. Let their right hand sides be F and G. 
We assume that G corresponds to a normalized quadratic vector field which 
is cr-distant from centers. This implies that F corresponds to a similar field 
which is at least 0.9cr-distant from centers. Let 



max IF - G\ < A, 
w 



max 
w 



dF 



dw 



where as before W = {\w\ < R} x S 1 , R = 0.01. Let e = e~ 2jTL R; clearly, 



m 
dw 
dz 



< %. Then the solutions wf and wq of the equations 



dw 

dz 



F and 



G with the same initial condition w(0) : \w(0)\ < | diverge on the 



segment < 9 < 2tt no more than 
(27) \w F {9) - w G (0)\ < 2TTAe 2wL . 

We apply (p7|) to our F and G. We have: 



A 



max 
W 



w\\ 



1 - wg x 



< 



mo 
96 



in W. On the other hand, L < 0.2 by (|19p . Hence, for any two solutions wf 
and wq with the initial condition w(0) and |io(0)| < e(0) = 0.0005, we have 

2ire 0An 

\w f (2tt) - w g (2tt)\ < — — — m < 0.4 m . 

9b 

Suppose now that w(0) = wf(0) corresponds to the solution wf for which 
\wf{2k) — wf(0)\ > 0.9 mo, and wc(0) = w F (0). Then \wg(2it) — wg(0)\ > 
and Lemma is proved. 



3. Upper estimate of the displacement of the Poincare map 

In this section we construct a neighborhood U\ of the set K\ where the 
Poincare map P\ of equation v\ is well defined. We give a lower estimate of 
the gap e between K\ and dU\, and find an upper estimate for fx = P\ — id. 
To this end, we find a universal gap between 5-tame limit cycles of quadratic 
vector fields that are K-distant from singular quadratic vector fields, and the 
curve 9 = 0. 
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3.1. The universal gap. A well known elementary property of quadratic 
vector fields ([I]) claims that any closed orbit of these fields that surrounds the 
singular point zero belongs to the domain 6 > 0. It is a simple consequence of 
the fact that any line has at most two contact points with a quadratic vector 
field. The boundary of this domain is given by the equation r = —l/g\{6). 

Lemma 12 (Second Main Lemma). No 5-tame limit cycle of a normal- 
ized vector field n-distant from singular quadratic vector fields intersects the 
curvilinear strip 



u p = {(e,r)eB x \ re 



1 ft 1 



for P 



5 u k 
10 10 ' 



9\{0) 9x{0). 

The proof of this lemma is technical. In the rest of this subsection we 
make the first step of the proof that makes the existence of the gap obvious. 
The estimates of the size of the gap are presented in Section |H 

Consider a zero isocline T: 

6 = 0, r ' 



gx(oy 

The restriction of f to this isocline equals 

r\r = o~ > H {v x ) = x i9x ~ fx- 

9x 

For the proof of Lemma [T2l we need a lower estimate of \H(v A)|rns A , 
First of all, we estimate from below the L2-norm ||^(va)|| 2 of H(v\) on 
S 1 = R/2vrZ. By JT2D, 

H(v x ) = Im Jih\. 

Note that H(v\) is linear with respect to v\. Let v \ = v s + u\ be the 
decomposition (JSj) for v\. For the singular vector field v s we have: H(v s ) = 0. 
Hence, 

Hivx) = H{u\) = Im fihx, 

where h\ = be~ %6 + ce~ 3l ° . 

Consider an arbitrary trigonometric polynomial H on R/27rZ. If H con- 
tains no complex conjugate monomials, that is, for any entry ae in ® + be~ md 
at least one coefficient is (i.e. ab = 0), then 

||Im H\\ 2 2 = ||Re H\\ 2 2 = -\\H\\ 2 2 . 

Indeed H = ^Za n e in6 implies that Re^ = \{Y,{a n e ine + a n e- n6 )), and 

consequently ||Reif|| 2 = ^S(l a «l + \®n\ ) = ^ll-^l^- ^he last conclusion 
holds because there are no cancelations in the sum for Re H , by assumption. 
The same argument proves the statement for Im H. 

Corollary 13. Forv\ which is n- distant from singular vector fields \ \H(v\)\\ 2 
H 
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Indeed, for equations, K-distant from singular ones, we have ||i?(fA)|| 2 = 

We got therefore a uniform lower bound for the Z^norm of the restriction 
r\r- It is now clear that a similar bound would exist for mmr\ rnB ^ X )- 
Indeed, zeros of r|r are located at the singular points of v \, and all the points 
of B(5,X) are at least <5-distant from these points. After mmr\j- nB ^x,S) is 
estimated, it is easy to prove that the lower boundary of the curvilinear 
strip Up has no contacts with the field v\. ^From this it follows that the 
5-tame limit cycles can not intersect ttr. The detailed proof of Lemma [T2l is 
completed in Section |H 



3.2. Construction of the larger domain U in the Growth-and-Zeros 
Theorem. Let 

S\ = S X \ £> e (A) 

and 



D = B(i - A>n i r£ ^^ 

For any A E A, consider a (/3<5)/32-neighborhood D' of the domain D C 
M + x S 1 in C x S . We will choose e in such a way that any orbit of v\ that 
starts in U £ x {0}, where U = U e is the e-neighborhood of S\, stays in D' 
while 9 ranges in [0, 2tc]. Let 



max 
D' 



dF x 



dw 



Then, by the Gronwall inequality, 



(28) E = ^e 

should be the desired one. Indeed, the largest <5-tame limit cycle keeps in 
D by Lemma [T2l Hence, all the orbits that start on S\ x {0}, keep in D 
by definition of S\. Then, for e from (|28j) . the orbits that start in U e x {0} 

would not quit D' for 9 E [0, 2tv]. Moreover, they will be close to the 

32 

real orbits starting at S\. Hence, the Poincare map for v\ is well defined in 
U e , and 

max I fxl = max \P\ — id I < 5~ l H . 

U e u 1 U e 1 1 ~ 32 

By Lemma the orbits that start in L^efA) s ^ a y i n Dr x S 1 as 9 ranges 
over [0, 2ir]. So, the set U\ = U £ U D 2e fx) * s a neighborhood of K\ in which 
the Poincare map of v\ is holomorphic, and 

(29) max |/ A | = M < 5~ x + 1. 
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3.3. The final estimate. We can now estimate the geometric factor in the 

Growth-and-Zeros Theorem . For this we need to get an upper bound for 

L, then a lower bound for e. 

To estimate L, we first get a lower estimate for the denominator in the 
dF 

relation (1181) for — — . We have: 
ow 

1 , & _ ox _! 



\w\d> <S +^« 26 



Now, estimate mhxjji \l + wg\\. If (w,9) € D' is such that \gx(&)\ < f > then 
\l + wg\\ > 1 — | • I > |. Suppose that |<?a(0)| is now greater than 4. Find 
a point (w',6) G D with |u/ — u;| < ||. Then 

. ^ . 1 , , . ^ J 86 ^ (58 

|1 + wg\\ > | ||sa| - |5a||w - w\ > (3- -4— > — . 

gx 4 32 8 

Hence,min£)/ |1 + wgx\ > 

Moreover, by LemmaEl Ai < 4<5 _1 . Hence, by (fl~8l) . 

L < 6145<5" 3 /?- 2 . 

We substitute this -L into (|28p and get the expression for e through 5 and (5. 
Note that the expression of (5 through 5, a, k is given in Lemma LT2l 
The intrinsic diameter D < 25~ l . Hence, 

— < 128r 2 /?- 1 e (1 ° 5 - 2)5 " 3/3 " 2 . 

£ 

2D 

This provides a double exponential estimate for the geometric factor exp . 

£ 

Note that for 5 < 0.1 and j3 < 0.1, increasing the factor 10 5 — 2 in the 
exponential by one will compensate well the division by the first factor. 
Finally, 

(30) ^ < ^105-1)5-3/3- 

We can now estimate the Bernstein index of fx- The numerator in (jlOf) is 
estimated in (|29p . The denominator is estimated in the First Main Lemma 
(Lemma [9]). We replace it by even smaller value: 

_26 

m = max \f\ \ > 10 s a. 

Kx ' ' ~ 

Finally, the Bernstein index of fx is: 

B Ux ,k x (fx) = log < log 2 - log 5 + ^ log 10 - log a. 

We see that this index, whose estimate took the main part of the work, is 
in a sense negligible in comparison with the geometric factor. Replacing of 
this index by | logcr| may be well compensated through the increasing by 1 
the exponential 10 5 — 1 in ([30]) . 



11 



YU ILYASHENKO AND JAUME LLIBRE 



Finally, by the Growth-and-Zeros Theorem we have: 

H(2,5,ct,k) < | log <r|e e 

Substituting here the value of (3 from Lemma H2] (which is not yet proved) , 
we obtain Theorem [5] 

4. Some lower bounds for trigonometric polynomials 
In this subsection we complete the proof of Lemma [12j 

4.1. Homogeneous polynomials of degree three. 

Lemma 14. Consider a real homogeneous trigonometric polynomial H of 
degree 3, that is, a homogeneous three-form on sin 9, cos 9 with real coeffi- 
cients. Let IRq, be the set of all real 9 that are at least a-distant from the 
complex rots of H . Then 

a 3 ,, „ 
min \H\ > — \\H \W. 

R a 24 

Proof. The polynomial H has three series of roots counted with multiplici- 
ties: 6j + irn, n £ Z, j = 1, 2, 3. Hence, for some real A, 

3 

H = Ajjsin(6»-0j). 

i 

Case 1. All 9j are real. Then 
(31) mm\H\>\A\(-) a 3 . 

On the other hand, 

The inequality: 2 2,5 /7r 3 ' 5 > 1/24 implies the lemma in Case 1. 

Case 2. One root 9\ is real, two others are complex: #2,3 = 9? ± iip, ij) 7^ 0. 
Then 

If = A\\shx{e-ej) = - s\n{9 - 9 1 ){ch2iP - cos 2(6» - ^)). 
1 

For any a 6 M, |6| < 7r, the following inequality holds: 
ch a — cos b> -a + — b . 



2 Vvr 

By assumption, ip 2 + (9 — (p) 2 > a 2 . Hence, once again we have ([3T1) . This 
proves the lemma in case 2. □ 
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4.2. Lower bounds for the distance to the roots. If two points of the 
disk r < 5~ l are at least <5-distant in Cartesian coordinates, then they are 
at least <5 2 -distant in the polar coordinates. If two points, one in the disk 
r < 5" 1 in M 2 , another in C 2 , are at least <5-distant in Cartesian coordinates, 
5 < 0.1, then they are at least |<5 2 -distant in complex polar coordinates. 

Proposition 15. Suppose that the point (0Q,r), r < S^ 1 and 9q £ [0, 2ir] is 
at least id 2 -distant from the singular points of the system (jlip with com- 



plexified r and 9 in the metric ds* 



(32) 

Then 

(33) 



r + 



9x(0 ) 



\dr\ 2 + \d6\ 2 , and 



5 2 



d(6 ,{H(v x ) = 0}) > 



5*_ 
100' 



Proof. By contraposition, assume that the converse to (|33|) is true. Then 
there exists 0\, zero of H(v\) such that 



h\ < a :- 



6*_ 
100' 



It may happen that 9\ is non-real. Take two extra points: b 
1 



1 



9x(9o) 



and c 



and let a = (0Q,r). Then, by ([32]) . 



Let L = max 
[60,61] 
Value Theorem 



$ 2 

a\ < — . 
1 ~ 2 



By assumption, \9q — 9%\ < a. Then, by the Mean 



\b-c\ < a\Jl? + 1. 



We now estimate L from above. Recall that g\ = Im h\, h\ = Ae lLp + 
Be-^ + Ce' 3 ^, \A\ < 1, \B\ < 2, \C\ < 1. Bv ([32]) and assumption r < S~ x , 
we have: 

1 <5 4 
|5a(^o)|>— ^ ' 



^2>5- — > 0.995. 

^ 



Now, by CED 

g x = Im h\, 



\g' x \ < \h' x \ < \A\e a + \B\e a + 3\C\e Sa < 7 



Then L < 85 2 . Hence, a\/ L 2 + 1 < — . Therefore, \a—c\ < la— 61 + 16— c| < 

6 



|^ 2 , a contradiction. 



□ 
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4.3. Proof of Lemma 1121 In order to prove that no limit cycle can cross 
11^, let us first check that the lower bound T~ of rig is a curve without 
contacts with the vector field (jlip . This lower bound has the form: 

(34) r-:r= - A (r,0) € B(X, S), 

Denote by S the minimal slope of the field (jlip on 

Ai + rh 



S 



mm 
r- 



dr 



mm 
r- 



On r we have: 

|Ai+r/ A | 



1 



(H( v \) - Pgxfx) 



1 + 



16/?). 



The points of -B(A, 5) are at least 5-distant from the singular points of system 
By the remark at the beginning of Subsection 14.21 points of T~ satisfy 
assumptions of Proposition [T5l Hence, for any 6 such that (r, 9) S T~ for 

some r, we have (|33p . Now, taking a = in Lemma [141 we conclude that 



100 



min \H(v\) 
r- 



> 



Q 

24 



fl2 



\H(v x ) 



10 6 • 24 



\H(v x ) 



By Corollary 1131 we get 



sl2 



min \ H(v\) 
r- 



> 



10 6 • 24\/2 



Hence 



Moreover, on T 



min I Ai + rf\\ 
r- 



> 



4/3. 



|1 + rg\\ 



-fax < 4/?. 



1 



At last, r|r > — . Hence 



S > 



1 

80/3 ~ 5' 

Denote by 7rr~ the projection of T~ to r 



along the r-axis; 7iT C 
{— g^ 1 < S^ 1 + /?}. We estimate the maximal slope of T~ . It is equal to 



max 

7iT- 



< 



6 



< 75~ 



mnvr- \g x \ 

The inequality S > s follows from the definition of (3 in Lemma [T2l 

We now prove that no <5-tame limit cycle that surrounds zero can cross 
n^. On the contrary, let a cycle 7 contain a point q G Hp. As 7 surrounds 0, 
it must enter and quit Hg. The connected component H q of that contains 
q is bounded by an arc 73^ of the curve (f34l) and by the part of dB(X,8). 
As S > s, the cycle can enter n 9 through 7^ (in positive or negative time) 
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but cannot quit U q through 7/3 l9 . Hence, it quits H q through dB(X, 5). This 
contradicts to the assumption that 7 is 5-tame and proves Lemma LT2l 
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6. The appendix 
In this appendix we provide the values of the a's, /3's and 7's of Lemma 

We shall compute the Poincare map P\ associated to the differential equa- 
tion (fl"4"|) in complex polar coordinates (u>, 9). Let P\ : {6 = 0} — > {8 = 0} be 
the Poincare map defined by the flow of system ()14|) ; i.e. P\ is the 27r-time 
Poincare map that brings an initial value of any solution r(#, x) of system 
(|14|) with initial condition r(0, x) = x on the half-axis {6 = 0} to the value 
of the same solution at 6 = 2ir, whenever defined. We know that the limit 
cycles surrounding the origin of system ([!]) correspond to real isolated zeros 
of the displacement function P\(x) — x. 

The power series expansion for the displacement function P\(x) — x as- 
sociated to a quadratic system (p} in a neighborhood of the origin is found 
in the following classical way. The right hand side of equation (I14j) may be 
decomposed in a power series in r with the ^-dependent coefficients: 

a 00 

(35) ■^r = E 1 *(*) w< > 

i=l 

where R\ = X\, 

(36) Ri(0) = (-mfx(9) - \igx(6)}g x (ey- 2 for i = 2, 3, . . . 

The modification of the Bautin result in [18J implies that the coefficients 
of the displacement map 

00 

(37) P x (x)-x = Y,a j (\)xi, 

i=i 

belong to the ideal generated by gj(X),j = 1, ...,4, where gj are the same as 
in Theorem [TJ For Ai = 0, the coefficients 07(A) are polynomial. 

We use the algorithm due to Bautin for computing explicitly P\(x) in 
powers of x up to order 7, see also [2]. We do the computations for the 
case Ai = 0; otherwise v?(9,\), which is necessary for computing vt(2tt, A) 
and consequently P\(x) in powers of x up to order 7, would need more than 
thousand pages. For doing these computations we have used the algebraic 
manipulator mathematica. 

We know that the series of (|35p converges if w is small enough, and that 
the solution w(9) of differential equation (|35l) satisfying the initial condition 
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w(0) = x can be expanded as 

oo 

(38) w(6,\)=Y,v i (0,\)x i , 

i=l 

where the Vi(9, A)'s satisfy the conditions 

(39) ui(0,A) = l and v;(0,A) = for i = 2,3,.... 

Substituting ([38]) in ([35]) . taking Ai = 0, and looking for the coefficients 
of the powers of x, we obtain the equations for determining all the v^s: 

^ = 
d9 

19 = 

dV3 

d9 



2v 1 v 2 R 2 + v\R 3 , 



= {2 Vl v 3 + vl)R 2 + 3vlv 2 R 3 + v\R A , 

^ = 2(viv 4 + v 2 v 3 )R2 + 3v l (viv 3 + vl)R 3 + 4vlv 2 R4 + vlR 5 

^ = (2v±v 5 + 2v 2 v/t + v 3 )R 2 + (3vfv 4 + §v 1 v 2 v 3 + v 2 )R 3 + 

2vl(2v lV3 + 3vl)R A + 5vfv 2 R 5 + v\R& , 

— = 2(viv 6 + v 2 v 5 + v 3 V4)R 2 + 
d9 

3(vfv 5 + 2viv 2 V4 + v\v1 + v%v 3 )R 3 + 



4^1 (V1V4 + 3?Ji?72W3 + V 2 )R4 + 

5wf (2u| + fif 3 )-R 5 + 6vlv 2 R 6 + ujl^ . 
All these differential equations are solved recursively computing an inte- 



gral with respect to 9 and taking into account the initial conditions (J39J). 
Thus, we get that vi(9,X) = 1, and 

v 2 (9, A) = -^(-3a 2 + 36 2 + c 2 ) + a 2 cos 9 - b 2 cos 9 - \-c 2 cos(3#)+ 

1 

a\ sin + 6i sin 9 H — ci sin(30) , 
3 

here we denote A = a\ +ia 2 , B = b\ + ib 2 and C = c\ + ic 2 . The expressions 
for v i(9, A) for i = 3, 4, 5, 6, 7 need approximately 1/2, 2, 7, 18 and 42 pages, 
respectively. Once we know Vi(9,X) for i = 3,4,5,6,7, evaluating Vi{2ir, A) 
we get the displacement function 

00 00 
(40) P x {x) = w(2n, A) = ]T Uj(27r, \)aP = ]T a, (A) aj> , 

i=i j=i 

with the explicit formulas for the polynomials aj(X),j = 1, ...,7. After that 
we decompose these polynomials in the ideal with generators gj, j = 1, 4. 
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This is done with the use of the manipulator mathematica again. The results 
of these computations presented below imply Lemma [101 The coeficients of 
the decompositions mentioned above are the following: 



ao = — 2-7T, 

fa = — -(9al - 96 2 a 2 - 6c 2 a 2 - 27vr6ia 2 + 276?. + 2lc 2 2 + 186?+ 
20c? + 66 2 c 2 - 276 2 7rai - 9ai&i), 

5"7T 

70= - T , 

71 = -.11 (300a 2 - 5586 2 a 2 - 240c 2 a 2 - 384?r6ia 2 + 528fo| + 204c 2 .- 

36a? + 2886? + 188c? + 1686 2 c 2 - 3846 2 vrai - I80161+ 
48aici + 24&ici), 

72 = __^_(2l60a| - 3606 2 a| - 1296c 2 a| - 25920vr6ia|+ 

171006|al + 27864cla| + 21600vr 2 6?a 2 . + 102606?a 2 .+ 
24648c?a 2 . + 72366 2 c 2 a^ - 259206 2 vraia| - 8280ai6ia 2 .+ 
345606 2 ?r6ia^ + 17280c 2 7T&ia 2 , - 4752aiCia 2 .- 
17406icia| - 342006|a 2 - 19824cia 2 - 34560vr6?a 2 - 
373686 2 c|a 2 + 46806 2 a?a 2 - 144c 2 a?a 2 - 334806 2 6?a 2 - 
9954c 2 6?a 2 + 172807rai6?a 2 - 384726 2 c?a 2 - 
20976c 2 c?a 2 - 38400vr6i c?a 2 - 228066^c 2 a 2 + 
345606^vraia 2 + 172 8 06 2 c 2 7raia 2 + 50406 2 ai6ia 2 + 
8280c 2 ai6io 2 + 432006 2 vr 2 ai6ia 2 - 604806|vr6ia 2 - 
41280c|7r6ia 2 - 172806 2 c 2 vr6ia 2 - 76086 2 aiCia 2 - 
1440c 2 aicia 2 - U628b 2 b 1 c 1 a 2 - 2\\2c 2 b 1 c 1 a 2 + 369006^+ 
13040(4 + 145806| + 11200cf + 161526 2 cf - 9720ai&?+ 
2640aicf + 1320&ic? + 707586 2 ,c 2 . - 41406?.a? + 648c?.a?+ 
21246 2 c 2 a? + 216006 2 .vr 2 a? + 507606 2 .6? + 42498c?.6?- 
1620a?6? - 255966 2 c 2 &? - 345606 2 vrai6? + 67074& 2 .c?+ 
24000cf c? + 120a?c? + 416706? c? + 152886 2 c 2 c?- 
384006 2 7raic? - 17400ai6ic? + 309966|c 2 - 604806^vrai- 
412806 2 c|vrai - YJ2mb\c 2 ^a x + 1080a?6i + 172806 2 7ra?6 r 
118806^ai6i - 19080c|ai&i + 5406 2 c 2 ai6i - 720a? c x - 
112206? a + 3690ai6?ci + 195 06 2 .aici + 2448c|aiCi+ 
13446 2 c 2 aici + 440286 2 6ici + 1224c^6iCi + 1500a?6ici+ 
13686 2 c 2 6ici). 
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